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ABSTRACT 

The evaluation, systematic analysis and numerical study of the semileptonic r-lepton de- 
cays with two mesons in the final state has been made in the frame of the standard model 
extended by right handed neutrinos. In the analysis, heavy-neutrino nondecoupling effects, 
finite quark masses, quark and meson mixings, finite widths of vector mesons, chiral sym- 
metry breakings in vector-meson-pseudoscalar-meson vertices and effective Higgs-boson- 
pseudoscalar-meson couplings have been included. Numerical estimates reveal that the 
decays r~ — > e~7r~7r + , t~ — > e~K~K + and r _ — > e~K°K° have branching ratios of the 
order of 10~ 6 , close to present-day experimental sensitivities. 



PACS number(s): 13.35.Dx, 14.60.St, 12.39.Fe 
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1 Introduction 



The neutrinoless r-lepton decays belong to the family of phenomena which, if ex- 
perimentaly confirmed, would unambiguously show that there exists physics beyond the 
standard model (SM). Specifically, the lepton sector would have to be modified. In the 
SM, these decays are forbidden, due to the fact the SM- neutrinos u e , and v T are exactly 
massless, the fact which follows from the doublet nature of neutrino and Higgs boson fields, 
left-handedness of the neutrinos, and chirality conservation. Neutrinoless r-lepton decays, 
if studied with sufficient accuracy, from the experimental point of view, are very promising 
due to the large momentum transfer involved In addition, the large mass of the r- 

lepton allows many decay channels. Therefore, SM (deviations from the SM) can be tested 
in a variety of ways. Experimental data on these decays constantly improve [0,0. The 
CLEO experiment H, has improved the previous upper bounds on 22 neutrinoless decay 
channels of the r-lepton by almost an order of magnitude. 

Neutrinoless r-lepton decays and many other lepton- number /flavour violating decays 
have been studied in a number of models, e.g. SU(2) x £7(1) theories with more than 
one Higgs doublet f|, leptoquark models ||, -R-parity violating supersymmetry scenarios 
0, superstring models with Eq symmetry ||, left-right symmetric models |J and theories 
containing heavy Dirac and/or Majorana neutrinos |T0| , pT| . Here, the models with heavy 
Dirac and/or Majorana neutrinos will be used to estimate the processes of interest. 

This paper is devoted to the analysis of semileptonic decays with two pseudoscalar 
mesons in the final state, denoted by r~ — > / =F P 1 P 2 - Together with papers |12| and [13[ 



it completes the analysis of the lepton number /flavour violating decays of the r-lepton 
reported by the CLEO collaboration H. In addition to the heavy- neutrino nondecou- 



pling effects [p^, p^Jl^Jl^JT6|] , finite quark mass contibutions, Cabbibo-Kobayashi-Maskawa 
(CKM) quark mixings and meson mixings already studied in the previous work [ 13| , this 
analysis includes vector-meson-pseudoscalar meson couplings, chiral symmetry-breaking 
effects, finite widths of the vector mesons and effective Higgs-pseudoscalar couplings. The 
hadronic matrix elements are derived in a few independent ways, in order to check the 
formalism used. 

For the evaluation of the leptonic part of the r - — > l T P\P2 matrix elements, the 
formalism and conventions of the model described in Ref. fl(| are adopted. The model is 
based on the SM group. Its neutrino sector is extended by the presence of a number (ur) 
of neutral isosinglets leading to heavy Majorana neutrinos (Nj). The quark sector of 
the model retains the SM structure. In couplings of charged and neutral current interac- 
tions, CKM-type matrices B and C appear HIU| , |l2" ,|rF|. These matrices satisfy a number 
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of identities, assuring the renormalisability of the model [0,0 and reducing the number 
of free parameters in the theory. These identities may be used to estabilish the relation 
between B and C matrices and neutrino masses, too. For example, in the model with two 



right-handed neutrinos, B and C matrices read [12] 
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where p 



rn 2 N2 /m 2 Nl , and s^ are heavy-light neutrino mixings |19| defined by 



nR 



(*2) 



El^l 2 - El- 8 ^ 



1.2) 



i=l 



The second equation (|1.2j) follows from the afore-mentioned relations for B and C ma- 
trices. In the theory with more than one isosinglet, the heavy- light neutrino mixing and 
light-neutrino masses {m Vl ) are not necessarily correlated through the traditional see-saw 
relation (s^) 2 oc m u Jm M . The (s^) 2 scales as (rrJ D (m]^) 2 m D )ii |]i7| , |T9|| , while light-neutrino 



masses depend on the matrix mcm^mj. If the condition m^m^m^ = is fullfiled, tree- 
level light-neutrino masses are equal zero, while (s^) 2 can assume large values. The light 
neutrinos receive nonzero values radiatively, but for reasonable tum values, their values are 
in agreement with the experimental upper bounds . Independence of the light-neutrino 



masses and the heavy-light neutrino mixings implies that (s^) 2 may be treated as free 



phenomenological parameters, which may be constrained by low energy data [19,20|. In 



this way, the following upper limits for the heavy-light neutrino mixings have been found 

2 < 1(T 8 . (1.3) 



K°L ) j K^lY < 0.015, K z L ) v u.wu, \ S L)\°L 

More recently, a global analysis of all available electroweak data accumulated at the CERN 



[s v £] 2 < 0.050, 



?(s? 



Large Electron Positron Collider (LEP) has yielded the more stringent limits |21 



< 0.0071, 



(s v £) 2 < 0.0014, 



< 0.033 (0.024 including LEP data) 



(1.4) 

at the 90% confidence level (CL). In this paper, the limits obtained in the Ref. |20) will 
be used because the results of the analysis in Ref. [2l|| depend to certain extent on the 



CL considered in the global analysis and on some model-dependent assumptions [12|. The 
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discussion on possible theoretical dependence of the upper limits, such as those in Eqs. ( |1.3| ) 
and (|1.4| ) , may be found in Ref . . 



The hadronic part of the amplitudes contains matrix elements of quark currents be- 
tween vacuum and a hadronic state. Vector and axial-vector quark currents are identi- 
fied with vector and pseudoscalar mesons through PCAC |E| and vector meson domi- 
nance |^,|4],[^| relations. The scalar quark current is expressed in terms of pseudoscalar 
mesons, identifying QCD and chiral-model Lagrangian. Intermediate vector mesons are 
described by the Breit-Wigner propagators with momentum- independent width p6|j27|j28|| . 
The vector-meson-pseudoscalar vertices are described by non-gauged U(3) l x U (3)r/U(3)v 
chiral Lagrangian containing hidden U(3)i oca i symmetry [25|, through which the vec- 



tor mesons are introduced. Both U{3)l x {7(3)fl/{7(3)y-symmetric and more realistic 
U{3)l x [7(3)i?/C/(3)y-broken Lagrangians |H| are used in the evaluation of the matrix 
elements. The gauge couplings of mesons are introduced indirectly through the quark 
gauge couplings in the above mentioned matrix elements of quark currents. 

This paper is organized as follows. In Section 2, the analytical expressions for branch- 
ing ratios of decay processes r~ — > e + PfP 2 ~ and r~ — > e~P 1 ~P 2 + /e P1P2 are derived. 
Technical details are relegated to the Appendices. Numerical results are presented in Sec- 
tion 3. Conclusions are given in Section 4. 



2 r- I'^PxPz 

In the model containing heavy Majorana neutrinos, there are two possible types of the 
semileptonic r-lepton decays into two pseudoscalar mesons 

1. T~ -> / /+ Pf P 2 ~ and 

2. T - -> I'-P^P? 2 , Qi + Q 2 = 0, 

where Pi and P 2 are pseudoscalar mesons, and Q\ and Q2 are their charges. Type (1) 
violates both lepton flavour and lepton number, and requires the exchange of Majorana 
neutrinos; henceforth these reactions will be refered to the Majorana type. Type (2) violates 
lepton flavour and proceeds via the exchange of Dirac or Majorana neutrinos; the appelation 
Dirac type will be attributed to these decays. Feynman diagrams pertinent to the Majorana- 
type and Dirac- type decays are given on Fig. 1(a) and Fig. 1(b), respectively. As mentioned 
in Introduction, only the decays with two-pseudoscalar final states, which are currently 
under experimental investigation, are considered. The decays with other two-meson final 
states could be calculated within the model, too, but they are phase-space suppressed, 
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they haven't been experimentally searched for, and they decay into the final states with 
more than two pseudocalar mesons. The complete calculation of such decays is much more 
involved than for the decays with two pseudoscalar mesons in the final state ||Z7[|. 

To start with, we consider the Majorana-type decays. At the lowest, fourth order in 
the weak interaction coupling constant, only tree diagrams contribute to the Majorana type 
decays. The chirality projection operators project out the mass terms of the numerators 
of the neutrino propagators. For that reason, only massive neutrinos contribute to the 
t~ — > / /+ P 1 _ P 2 ~ amplitude. Since the H^-boson and heavy neutrino masses [ 10 1 are much 
larger than the energy scale at which quarks hadronize to mesons, their propagators may 
be shrunk to points so as to form an effective amplitude depending only on one space-time 
coordinate: 

S(r- - V+P X P 2 ) = £ V* d V* d £ - 7 5 K / *xe~«*-rt* 

2M w o,6=l i=i mjv * J 

{P^P^\d a {x)^{\ - lh )u{x)d b {x) llx {l - -y 5 )u(x)\0), (2.1) 

where aw = a em / sin 2 9w ~ 0.0323 is the weak fine-structure constant, Mw is IV-boson 
mass, V u d a are CKM matrix elements, are heavy neutrino masses and u(x) and d a {x) 
are quark fields for u, d and s quark (di = d and d 2 = s). A more reliable calculation would 
also include the QCD corrections of four quark operators in eq. fl2~T| ) (they introduce new 
quark operators, and mixing of all quark operators), along with a renormalization-group 
analysis of their coefficients [plj^2| . Since such refinements will not alter our conclusions 
concerning the magnitude of the amplitude, they will be ignored. 

The hadronic matrix element may be evaluated using vacuum saturation approxi- 



mation and PCAC Vacuum saturation approximation [p2| , |33| allows to split the matrix 
elements involving four-quark operators into matrix elements of two-quark operators. The 
two-quark operators forming axial-vector currents may be combined into the currents hav- 
ing the same quark content as the produced pseudoscalar mesons, P, A^(x). The matrix 
elements of the currents A^(x) are evaluated using the PCAC relation f2"2| : 

(0\A?(x)\P'(p P ,)) = 6pp.V2fp.rfe-*'", (2.2) 

where f P i is the decay constant of pseudoscalar meson P'. The Kronecker symbol Spp. 
assures that the matrix elements ( p.2|) give the nonzero result only if the final state quantum 
numbers match those of the axial- vector current. Following the above procedure, one 
obtains the expression for the generic matrix element of r~ — > l' + P{ P^ process 

tIt" >. /'+ p~ p~\ — i8a w iT 2 fpxfpt D * p* _J 

6 M W i=l m N l 

x (PpiPp 2 )%(! - lb)u T . (2.3) 
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The corresponding branching ratio reads 



v b b Mw 2 r 



l\2 



i>(m—m') 

dt u, 

'(mi+m2) 2 



36r T m 3 M^ 

(2.4) 

where S is the statistical factor, equal to 1/2 if two equal pseudoscalars appear in the final 
state, and uj is a phase-space integral of the Mandelstam-variables dependent part of the 
square of the amplitude which is defined in Appendix C. 

Now we turn to the Dirac-type decays. The scattering matrix element of r _ — > l'~ P\P 2 
receives contributions from 7-exchange graphs, Z-boson-exchange graphs, box graphs, 
Higgs-boson(if )-exchange graphs and W /+ -boson-lV~-boson-exchange graphs, 

S(r- _> l'- Pl p 2 ) = 5 7 (r- -> 1'-P 1 P 2 ) + Sz(t~ -> l'~ PiP 2 ) + S Box (r- -> l'~ P^) 

+S H (r- - rPxP 2 ) + S w - w+ (r- - rPxP 2 ). (2.5) 

The 7, Z-boson and Higgs-boson amplitudes factorize into leptonic vertex corrections and 
hadronic pieces. The loop integrations are straightforward. The hadronic parts of the 7- 
and Z-boson amplitudes consist of the vacuum-to- vector-meson matrix element of the local 
vector and axial-vector quark current (only vector quark currents have nonzero contribu- 
tions, since only vector mesons decay into the two-pseudoscalar-meson state), a propagator 
of the vector meson and the vector-meson- P\-P 2 vertex. The hadronic part of the H 
amplitude contains vacuum-to- P\-P 2 matrix element of the local scalar quark current. Ex- 
ploiting translation invariance, the phases that describe the motion of the meson(s) formed 
in a vacuum-to-hadron matrix element may be isolated. Therefore, only the space-time 
independent hadronic matrix elements remain. These phases assure four-momentum con- 
servation. The 7, Z-boson and Higgs-boson amplitudes read 

S y (r~ ^1'-P 1 P 2 ) = (27t)U^ (p-p'- Pl -p 2 ) 

E^"(r - l'V°) iS V04 M T"(V° - P X P 2 ) 

S z (t- -+l'-P 1 P 2 ) = (2n) 4 5^(p-p'-p 1 -p 2 ) 

£T£(r - l'V°) iSyo^(q) T»(V° -> P\P 2 ) 

yo 

S H (r- ^1'-P 1 P 2 ) = (2n)*5^(p-p'-p 1 -p 2 )T H (T^l'P 1 P 2 ), (2.6) 

where p,p',Pi and p 2 are the four-momenta of r, I', Pi and P 2 , respectively, J2v° * s a sum 
over vector mesons that appear simultaneously in T^ z and T U (V° — > P\P 2 ) amplitudes, 
S v o !flu (q) is a constant-width Breit-Wigner propagator |26]j27|j28|| of the vector meson Vq, 

— n -4- q ^ v 

y^u 1 ~mT~ 

Sv °^ M) = g 2 - M? + i M v0 T v0 ' (2 ' 7) 
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T U (V° — > P1P2) multiplied by the polarisation vector, (q), gives a V°-Pi-P 2 vertex, 
which may be read from the Lagrangians ( A.l ) and ( A . 1 1|) , T^ z (r — > i'V" ) are 7- and im- 
parts of the T-matrix elements for the r — > Z'^/ reaction [T^|, from which a polarization 
vector of the l/ meson is removed, 
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P|< u,, 7 "(l - 75 )u T 



x((y |w(0)7 A1 (l-75--<)«(0)|0) 



-(y°M(0) T ,(l- 75-3^)«i0) 
-(V°|.(0)7,(l-75-^)K0)|0)Y 



;2.9) 



and Th(t — > VP1P2) is the T-matrix element of the r — > P1P2 reaction, 



T H (r^l'P 1 P 2 



w 



{mu v (l + 7 5 )m t P/ + m'«i»(l - 75 K Gjf) 



(P 1 P 2 |m„M(0)M(0) +m d J(0)cZ(0) + m s s(0)s(0)|0) 



(2-10) 



In Eqs. ( f2T7] - |2.iO| ) m T , m', Mh, rn u , ma and m s are masses of the r, /', Higgs boson, u, d 
and s quarks respectively, sw = sin 6 W is sine of the Weinberg angle, and represent 
r — > Z'7 and r — > I'Z loop functions, respectively, multiplied by corresponding gauge-boson 
propagators, j^ m (0) is quark electromagnetic current, and Kf(0) and A^(0) are vector and 
axial- vector quark currents for quark-Z-bozon interaction. The loop form factors FJj and 
Gjj may be found in Appendix B and FT. 1 ' and P| z in Eq. (|2~T 



in Ref. 12 



The box and W + -W~ diagrams are more involved as they contain bilocal hadron 
currents. In the case of the box diagram, the bilocality problem can be overwhelmed since 
the two VF-bosons in the loop assure the high virtualities of the loop momenta. That 
allows one to approximate the loop-quark propagator with the free quark propagator, and 



to replace the bilocal vector and axial- vector current operators with the local ones [13|. As 
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in 7- and Z- amplitudes, only the vector quark current operators contribute, giving rise to 
the vector mesons, which decay into the two-pseudoscalar-meson final state. In this way 
one arrives at the following expression for the box S'-matrix element 

S Box (r- - l'-PxP2) = (2n)' l 6^\p-p' - Pl - p 2 ) 

HTLAi - l'V°) iSyo^(<l) T V {V° - P 1 P 2 ), (2.11) 

yo 

where T Box (l — > l'V°) is the box part of the T-matrix element for the process I — > l'V° [12 
from which the polarization vector of the vector meson, V°, is removed, 

T Box (l^l'V°) = n ox {l ^ l'V°)el\q) = ^ )M (f°|y/ o ^(0)-^^(0)|0) 
" E ^„ A *<t^|^«'***(0)-Af»A*(o)|o> 



d„h=d,s 



1 a 



2 

w - 



16M^ 



un^l ~ 7sX F T lT (V°|m(0)7"(1 - 7 5 )«(0)|0> 



E n™ adb (V°\d a (0mi- 75)4(0)|0)1, (2.12) 



=d,s 



where L Boxm i are box loop functions, and V^ ox ' qq '(0) and A^ ox,qq ' (0) are the corresponding 
vector and axial-vector quark currents in an r — ► I'qq' amplitude. The loop form factors 
F T B V d x adb and F T B v ox u are defined in Ref. |3| . 



As in the r~ — > l~ P{ P 2 ~ amplitude, the W-bosons in the W /+ -W / ~-exchange diagram 
may be shrunk to points. So, an effective amplitude depending on two space coordinates 
is formed. The chiral projection operators extract the momentum dependent parts of the 
numerators of the neutrino propagators, so that both heavy and light neutrinos contribute. 
The heavy-neutrino propagators could also be shrunk to a point, and, therefore, the cor- 
responding amplitudes depend on one space-time coordinate. By contrast, light-neutrino 
contibutions cannot be reduced from the bilocal to a local form. To enable the comparison 
of contributions of heavy and light neutrinos, all contributions to the transition matrix 
element are written in their bilocal form, 

dH 



S{r-^l'- Pl P 2 )= ^ E V: d V udb ^B mi B: Ni jd^y- ■ 

W d a j,=d,s i=l v / 

r m N 

x (P 1 P 2 \ u{y) lv {l -K)d b (y)d a (x)^(l -^)u{x) |0). (2.13) 



As I 2 < m 2 . and the lightest heavy- neutrino mass exceeds 100 GeV [[L(|, the local (heavy- 
neutrino) terms are supressed at least by factor 10~ 4 relatively to the nonlocal (light- 
neutrino) terms. Therefore, one can safely neglect them. 



S 



The amplitudes ( |2.6| ), ( |2.11| ) and ( |2.13| ) comprise three types of hadronic matrix 
elements: (V°\q(0)^q(0)\0) , (P 1 P 2 \m g q(0)q(0)\0) and {P x P 2 \u{x) llM d a {x)d b {y) lv u{y)\Q). 

The evaluation of the (V r °|?(0)7 A1 g(0)|0) matrix element proceeds as follows. The two- 
quark operator (?(0)7 M g(0) is expressed in terms of vector currents, V^, having the same 
quark content as the produced vector mesons, V°. Exploiting the vector-meson dominance 
relation |23|], correlating a vector-meson field V^(x) and vector current, V^, having the same 
quark content as V^(x), 



2 



one arrives at the expression 
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1^ (X)\V (P V o)) ^ Sy l y - 7 =^—£y 0|x (pyO,Xy )e 



■^P i rO^ 



(2.14) 



(2.15) 



The Kronecker symbol S v , v0 , assures that the matrix elements give non-zero contributions 
only if the vector-meson quantum numbers match those of the vector current. 

The (PiP 2 \ J2q= u ,d,s m q(l(^) ( l(^)\^) matrix elements may be evaluated comparing the 
quark sector of the SM Lagrangian, and the corresponding effective chiral Lagrangian, 
contained in the first and second curly bracket of Eq. ( |A.1| ), one obtains the expression for 



the scalar two-quark current in terms of pseudoscalar fields [p3 

1 



q(x) l q(x)° = ~-jlr 



U(x) + U(x) ] 



>j 



(2.16) 



where U(x) = exp(2iTi(x) / / OT 



7T (X) 



T a Ti a (x), 7r a (x) are pseudoscalar meson fields, T a 



X a /2, X a are the Gell-Mann matrices and 

2m! 



2m 2 K0 



2m 2 K , 



m d + m u m d + m s m u + m s 



(2.17) 



Exploiting Eq. ( |2.16| ), one can write the H — q — q part of the Yukawa Lagrangian in terms 
of pseudoscalar fields 



C 



Hqq 



9w 



2M W 
9w 



H(x) m q q(x)q(x) 



q=u,d,s 



H(x) 



AM W 

+m 2 K0 K°{x)K°{x) + 



ml(ix-{x)Ti + {x) + 7r°(a;)7r (a;)) + m 2 K+ K + {x)R-{x) 
2^2, 



■{2ml ~ m K+ ~ m^ )r]i(x)?7 8 (x) 



+\( m K+ + m K° + m l)v!(x) + \{2m 2 K+ + 2m 2 K0 - m^)^(x) 



(2.18) 



9 



where H(x) is the Higgs field and ir~(x), tt + (x), tt°(x) etc. are pseudoscalar-meson fields. 
Replacing the fields r) s (x) and T]i(x) by physical fields r](x) and rj [x) given in Table II, one 
obtains the set of i7-boson-pseudoscalar-meson couplings. 

The evaluation of the {PiP 2 \u(x)'j^d a (x)db(y) r yuu(y)\0) matrix element is, in its full 
complexity, a highly nonpertubative problem due to the nonlocality of the four-quark op- 
erators. The one-loop pertubative QCD analysis of the W + W~ diagram shows that the 
corresponding amplitude has strong IR divergencies, but no UV divergencies, even if W- 
propagators are shrunk to points. That suggest the evaluation of the matrix element in the 
model which is valid at very low energies, the gauged U(3)l x U(3)r/U(3)v chiral model 
with pseudoscalar mesons coupled to the WS gauge bosons. The calculations in the chiral 
model show that the contributions to the amplitude come only from the diagrams with 
pseudoscalar mesons emitted from different space time points. In the quark picture that 



would correspond to splitting of the hadronic matrix element ( |2.13|) into two vacuum to 
pseudoscalar-meson matrix elements of the two quark operators, 

{PiP 2 \u{x)^^{\ - ^)d a {x)d h {y)^ u {\ - j 5 )u(y)\0) 

« (P 1 \u(x)j^ 5 d a (x)\0){P 2 \d b (y)^ l5 u(y)\0) + (Pi <-> P 2 ) 

= 2f P J P2 5 PlP{udc Jp 2 p {dbU c ) e^ x e^yp 1 ,p 2u + (P 1 <-> P 2 ), (2.19) 

where P(ud^) and P(dt,u c ) are pseudoscalar mesons having quantum numbers of the com- 
binations of quarks in brackets (q c is symbol for antiquark). Both chiral model approach 
and quark model approach, in which Q2.19|) is assumed, give the same results. Although the 
obtained result is appealing, one must have in mind that chiral models work for momentum 
transfers ~ 1 GeV 2 . Therefore, it is worth to compare this results with results obtained by 
some other method, e.g. sum rules. In the sum rule approach, it is quite unlikely that one 
can split the matrix element as in Eq. (|2.19|) , and consequently the quarks coming from 
the different space-time points are expected to form the (neutral) pseudoscalar mesons, 
also. That somewhat lessens the value of the approximation (|2.19|) . Unfortunately, the 
matrix element with two light pseudoscalar mesons in the final state cannot be treated by 
usual sum-rule techniques as in case of processes with only one light pseudoscalar meson in 
the final state, as for instance in D* — > Dtt decays [Sl|, because of complications of large 



distance strong interactions. The approximation (|2.19|) will be used here, because from 
phenomenology it is known that such approximation can hardly fail the correct value of 
the amplitude by a factor larger than 5, and because chiral model calculation suggest that 
approximation. 
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Following the procedure outlined above, one obtains the expression for the generic 
T(t~ — > / / ~P 1 P 2 ) matrix element 



T(t~ - 1'-P 1 P 2 ) = u n ,(l - 7 sK «' Pa (Pi - P2T + B% P2 qn 

+uv^f- (m(l + 75) + m'(l - 7s )K C7g' ft (r - p 2 f 
+u v {\ + 7 5 )w T £>g' ft + tZ,/(l - 7 5 )w T £ P fp 2 

&(^~ A) A(l - J5)UtF£ P2 (2.20) 

The first two terms belong to the 7, Z-boson and box amplitude, the third and fourth to 
the Higgs-boson amplitude, and the last one to the W + — W~ amplitude. The composite 
form factors A P { P2 , B T P l [ P2 , Cpfp 2 , D T P l [ P2 , E T P \ P2 and Fpf Pa read 

A p\p 2 = -T,Plw((l)Cv"p 1 p 2 i(a^ +b T ] l'o) 
v° 

2 2 

itV Y^„V° f„\n ;f„Tl' , ltI' \ m l ~ m 2 



B p\p 2 = J2PBw(<l) C! VOp 1 P 2 K a VO+ b 



y0) Ml 



V° Vo 

Cp[ p 2 = X! P V bw ( l) Cv Pi p 2 i c v° 
v° 

ryrV _ la W M HP 1 P 2 t?tV 

v _ iot^ Mjj 0Pi p 2 v 

~ ISM*, Ml mGn 
2 2 

• a H' f T/ T/* f f T?tV 



F PiP 2 - i A #4 V nd a V* d Jp 1 fp 2 F^ v - w+ , (2.21) 
1V1 W 



where 



Pew = T 5—^ 7T- (2-22) 

is a denominator-part of Breit-Wigner propagator for a vector meson V° (|2.7|) . Cyop 1 p 2 are 
— Pi — P 2 couplings defined by Lagrangian (|A.1|) , Oy , 6y ,and Cy' a are composite form 
factors for r — > Z'^ decays found in Ref. [0 and listed in Appendix B, and F^- w + is the 
tree-level form factor, 

F w-w+ = 7 ^^Bl'NiBrNi- (2.23) 

Here few comments are in order. 



1. From the structure of the total amplitude ( |2.20| ), one can easily find 
which of the amplitudes T 7 , Tz, Tp x, Th and T w - W + give the dom- 
inant contribution. The amplitudes T 7 , Tz and T Pox contain a com- 
mon factor (ia^ v /16M^ v )(g p7T7r /'y v ). In place of that factor, in the 
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amplitudes T# and T w - W + are factors (ia 2 v /16M^ v )(M HPi p 2 /M H ) and 
{ia%?K 2 /M^){f P JpjMlr)V ud J* db Y l N i BvN i B* N ., respectively. The amplitudes 
Tz and Th contain loop form factors behaving as the square of the heavy neu- 
trino mass, m 2 N , in the large-m^ limit, T 7 and Tbox have In asymptotics in 
that limit, and T w - W + is almost independent on m N . Approximating roughly all 
momenta of outer particles with r-lepton mass, one obtains approximate ratio 
of the magnitudes of the amplitudes 

Ml 



r '/■ . r QpTYTV tptV . IV1 HPiP 2 F Tf 
J--y,Z,Box ■ J-H ■ J-W-W+ ~ P z . —2 P H 

lp M H 



l^ 2 ^J^V uda V: db £ B VNi B* Ni (2.24) 



» ,f 2 ua b 
1V1 W Ni 

For heavy-light neutrino mixings {s v £) 2 = 0.01, (s^ 1 ) 2 = and {s u £) 2 = 0.05, 
Fp' and F H V assume values —0.01 and 0.01, respectively, for = 100 GeV, 
and values —1.6 and 2.2, respectively, for maximal value of allowed by the 
pertubative unitarity relation [see Eq. (|3.1| ) below], tun = 3700 GeV. Putting 
these values into Eq. ( |2.24| ), one finds that the T w - W + and T H amplitudes are 
six to four and four orders of magnitude smaller than the T 7 ^ z ,Box amplitude, 
respectively. The numerical study of relative T lt z,Box, Th and T w - W + contri- 
butions to the r~ — > V P1P2 branching ratios shows that the T J} z,Box amplitude 
participates even more than forseen by this rough estimate. Therefore, one can 
safely neglect H and W~W + contributions in the expressions for the largest 
branching ratios. Since within approximation ( |2.19|) only Tjj amplitude partic- 
ipates to t~ — > l'^7T 7i /l'~i]ri/l'~riri' channels, it will be kept for illustration of 
magnitudes of corresponding branching ratios in Fig. 2. 

2. As mentioned in Introduction, the hadronic matrix elements are evaluated using 
the non-gauged U{3)l X U(3)n/U(3)y Lagrangian containing hidden U(3)i oca i 
local symmetries. The effective gauge-boson-meson couplings are introduced 
through the gauge-boson-quark couplings and PCAC (|2.2|) and vector meson 
dominance ( |2.14| ) relations. The corresponding effective Lagrangians for vector- 



boson-7 and vector-boson-Z interactions read 



Ac w 



9W ry^p Q , m l( C V C 2W Sy 



. ml (SyC 2 W Cy \ 



(2.25) 
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where sy = sin6V and cy = cos9y. The 7, Z and W~W + amplitudes could 
be also evaluated using the gauged version of the U(3)l X U(3)r/U(3) v chiral 
Lagrangian with hidden U(3)i oca i symmetry. Both approaches give the same 
results for these amplitudes. That follows from the comparison of the effective 
Lagrangians ( |2.25|) and the corresponding terms in the gauged chiral Lagrangian 
flOP- Identifying 

agf l = £ = < (2. 26) 

27 P 27^ 27^ 

the Lagrangians ( |2.25| ) and the coresponding parts of the Lagrangian (|A.1|) be- 
come equal. This identification is justified numerically. The same type of iden- 
tification for jy-boson-pseudoscalar-meson couplings is trivial, because both 
approaches use the same hadronic parameters, pseudoscalar-meson decay con- 
stants. The indirect way to evaluate hadronic part of the amplitudes was chosen 
because the Tbox and Th amplitudes do not have their chiral model counter- 
parts. Moreover, this approach enables one to use the experimental values for 
the meson masses and branching ratios. In the chiral model, they are determined 
by the symmetries of the model. 

3. The chiral nonlinear Lagrangian based on the U(3)l x U(3)r/U(3)v symme- 



try (without hidden symmetries) describes well the treshold processes [ [2"8] , f2~9[j 
with pseudoscalar mesons in the final state only, i.e. amplitudes of vanishing 
pseudoscalar momenta. To comprise the dominant two-pseudoscalar channels 
of the final state interactions which swich on at higher energies, vector mesons 
are introduced . One of the most common ways to include the effects of pres- 
ence of vector mesons into the low energy chiral model amplitudes is to multiply 
them with the Breit-Wigner propagators normalized to unity at zero-momentum 
transfer. The constant-width normalized Breit-Wigner propagator has the fol- 
lowing form [p^,p7|,pS 



Ml - zM f/ r f/ 

Y -^—, (2.27) 

Ml-t-iMyTy' 1 J 

where M v and T v are vector meson mass and decay width, respectively. The 7, 
Z and box amplitudes obtained in the formalism of this paper have almost the 
same structure, 

T^ Box = L^^iiP^yolV^*™ - Al' z > B °*\Q) x K^ Box 

M 2 ~ 

x_ ^YH (2.28) 

M v0 -t- iM v oT v o K 1 



13 



where L^ ZBox are 1°°P P ar ts of the T — * I' P\Pi amplitude defined in Eqs. (|2T 
I| and gig ), K r 

z,Box are factors containing coupling constants (K^ = —ie, 
K z = -ig w /Ac w and K Box = 1), and ((PiP 2 )yo I ^7' z ' Bos - AJ z ' Soa; |0) comprise 
products of a vacuum-to-vector meson amplitudes of a quark current devided by 
square of the vector meson mass, a denominator of the vector-meson propagators 
and a vector- meson-pseudoscalar- meson vertex. The factor M v , which devides 
the vacuum-to- vector meson amplitude of the quark current, is extracted from 
the composite form factors for the r~ — > l'~V°, Oy , by' and Cy' , and is assigned 
to the vector meson propagator. The low energy limit of the matrix elements 
(PiP2\V / J' z,Box — AJ l ' z,Box \0) may be derived from the kinetic part of the chiral 
part of the Lagrangian ( |A.1|) , (f%/ '4)Tr(c^L? d^W) , identifying the quark vec- 
tor currents with the corresponding pseudoscalar-meson vector currents which 
may be found in Appendix A. These low energy limit amplitudes coincide with 
the corresponding amplitudes in Eq. Q2.28 ) for zero momentum transfer if the 
replacement 

M' v — My - iMyYy (2.29) 

is made, if 

1p = lw = 14> ( 2 - 30 ) 

^ = 1 (2.31) 

is made. The equality of the factors 7y is a consequence of the U(3)i x 
U(3)r/U(3)v symmetry, and relation ( |2.31|) is nothing but the famous 
Kawarabayashi-Suzuki-Riazuddin-Fayazuddin relation |^6j. Therefore, only the 
replacement ( |2.29| ) has no natural explanation. It will be included "by hand", 
by replacing 

My Ml - iUyYy 

v v v v ^ 2 _3^ 



and if the identification 



V2 lv V2 lv 
in the vector- meson-dominance relation fl2.14|) . 

4. The Lagrangian f |OD has U(3) L x U(3) R /U(3) V symmetry. The breaking of that 
symmetry will be introduced in the way of Bando, Kugo and Yamawaki 30] by 
adding extra terms in the Lagrangian [compare Eqs. ( |A.1| ) and (|A.11|) 1 and by 
renormalizing the pseudoscalar fields. In that way, the hidden U(3)i oca i symme- 
try, which becomes dependent on U(3)l x U(3)r symmetry through the gauge 
fixing, is also broken. Since the Bando et al. Lagrangian is not hermitean, the 
Lagrangian in Eq. (|A.11|) is written as half of the sum of their Lagrangian and 
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its hermitean conjugate. Assuming the ideal mixing between S77(3)-octet and 
S"?7(3)-singlet vector meson states, Oy = arctan(l/A/2), Bando, Kugo and Ya- 
mawaki obtained the following relations between pseudoscalar decay constants, 
vector meson masses and vector meson gauge coupling constants 

Sk 



'1 2 z r- 

m P = m u = a 9 L 
m 2 p ml y/2ml g 



2 2 2 f2 " b K* uu^ 

p m u - ag In - i+ - (i . Cy Y. 



(2.33) 



where Ca and Cy are breaking parameters appearing in the Lagrangian ( A. 11 ), 
and g yp , g 1(jj and g^ are gauge-boson-vector meson coupling constants which 
may be read from the Lagrangians ( |A.1|) and ( [A.ll|) . Replacing the expressions 
for the gauge coupling constants from Eq. ( |2.33| ) with the corresponding expres- 
sions in the Lagrangians ( 2.25|) into the third of Eqs. (|2.33|) , one obtains again 



Eq. ( |2.30|) . Therefore, if the ideal mixing between 5'?7(3)-octet and S'L r (3)-singlet 
vector mesons is assumed, the equality of 7yo-s is preserved after the symme- 
try breaking. In this paper, the ideal mixing condition is relaxed: the mixing 
angle 9y is evaluated from the experimental meson masses using the quadratic 
Gell-Mann-Okubo mass formula. 

Keeping in mind the above comments, one can derive the corresponding expression 
for the branching ratios from the expression for the generic r~ — > l'~ PiP 2 amplitude: 

I r(rn—m') 2 rsi~ 

„n2 



/-(m— m') 
(mi+m 2 ) 2 

rl' |2 st atV riTl'* i l „ \ |2 



/ dt [a\A P [ P f + P(A P [ P2 B P [ P2 + h.c.) 



+l\B P [p 2 \ 2 - SiA^C^ + h.c.) - e\C P [ P2 \ 
H{A P U<k + ^E P %)+h.c.] 

+v{B P U<k + ^E P %) + h.c.] 
+$(C P UD P % + ^E P %) + h. c : / 
M\D t p[p 2 \ 2 + KpJ 2 ) + ><<pXp 2 + h.c.)}, (2.34) 



where integration boundaries sf and parts of the square of the amplitude depending on 
the momentum transfer variable t, a, (3, 7, 5, e, £, m 1 and k may be found in Appendix 
C. 
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3 Numerical results 



In the numerical analysis, the extension of the SM with two heavy neutrinos is assumed. 
The description of the model and the relevant formulas for B and C matrices may be found 
in Introduction. The additional parameters of the model are three heavy-light mixings, 
s^, and two heavy neutrino masses, m Nl and m N2 . The upper limits ( |1.3| ) and ( |1.4j ) 



experimentally constrain the mixings s^, while the upper bound on heavy neutrino masses, 



[£(^)T\ P > 1 (3-1) 



may be obtained from the perturbative unitarity relations |I2| , |I5[ft7|| . The experimental 
upper bound limits (|1.3| ) suggest that either s v £ or s v £ is approximately equal zero. Here 
will be assumed that s v £ ~ 0, and, therefore, only t~ — > e T PiP 2 decays are considered. 
The results obtained for s u £ ~ case, that is for r~ — > p^P\P2 decays, almost coincide 
with corresponding s£ ~ results, and it is superfluous to discuss them seperately. The 
t~ — f e T P\P2 decays depend on new parameters of the model, s v £ y s v £ and m^., as well 
as on a whole set of quark-level parameters and meson observables: CKM mixing angles, 
quark masses, mixing angle between octet and singlet vector-meson states, meson masses 
and decay widths, pseudoscalar-meson decay constants, constants describing the coupling 
strength of vector mesons to the gauge bosons and vector-meson-psudoscalar-meson cou- 
pling constants. In calculations, the average of the experimental upper and lower values 
for CKM matrix elements are used, and the quark masses 

m u = 0.005 GeV, m d = O.OlOGeV, m s = 0.199 GeV, 

m c = 1.35 GeV, m b = 4.3 GeV, m t = 176 GeV, (3.2) 

cited in Refs. fl38| , |3~9"|l . The masses off all quarks are kept in evaluation of matrix ele- 
ments, since t and c quarks give comparable contributions to some amplitudes. The mixing 
angle between singlet and octet vector-meson states is not taken to be equal to the ideal- 
mixing value, 9y = arctan(l/-y/2), but is either determined from the quadratic Gell-Mann- 
Okubo mass formula, or treated as a free parameter. For pseudoscalar decay constants 
f n ± and /i<-±, appearing only in the W + W~ amplitudes of r~ — > e+P^P^ decays and 
t~ — > e + P{ P2 amplitudes, the experimental values are used 

f n ± = 92.4 MeV, f K ± = 113 MeV. (3.3) 

The constants jy, describing the coupling strengths of vector mesons to the gauge bosons, 
are either extracted from V —>■ e + e~ decay rates 

7p o = 2.519, 7w = 2.841, 1(j> = 3.037, (3.4) 
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or estimated using SU (3)-octet symmetry: 7^0* = 7 p o. Notice that the equality of 7yo-s 
predicted by U{3)l x U(3)r/U(3)v symmetric chiral model and by U(3)l x U(3)r/U(3)v 
broken chiral model is reasonably satisfied. The decay rates of vector mesons, involved 
through the vector-meson propagators, are taken to be equal to their experimental total- 



decay-rate values [pq| , and are not treated as momentum dependent quantities [p7[| . The 
p-TT-ir coupling is derived from the p — ► 2tt decay width, while the other vector-meson- 
pseudoscalar-meson couplings are fixed by one of the chiral models described in Appendix 
A. It is visible from the above that, whenever possible, the parameters were extracted from 
experiment and model dependent relations determining them were relaxed. 

In this paper, 17 r~ — > e T P\P2 decays are studied numerically. For orientation of the 
reader, decay widths of all 17 reactions are plotted in Fig. 2 as functions of = \rriN 2 
for upper bound values of heavy-light neutrino mixings ( |1.3| ). Concerning the de- 
pendence, the decays can be split into four groups: — > e~7r + 7r~ /e~K + K~ /e~K°K°, 
T~ -> e-7r + K-/e-n-K + /e-TT K /e-7r K /e~K r]/e'K r]/e~K°r]'/e-K 7]', r~ -+ 
er 7r° '7T /e~i]ri/e~riri' and r _ — ► e + 7T~7T~ / e + 7T~ K~ / e + K~ K~ . Only the decays of the first 
group are interesting from the experimental point of view and receive contributions from 
all five r~ — > e~PiP 2 amplitudes [see Eq. (|2.5|) ]. The others are suppressed by at least 
8 orders of magnitude relative to the first group of decays for various reasons. The mem- 
bers of the second group are Cabbibo suppressed, and only box and W + W~ diagrams 
contribute to them. The decays of the third group originate from the //-amplitude and 
are suppressed by the factor {M\ PxP J M^) 2 from Eq. ( |2.24| ). The last group belongs to 
the Majorana-type decays, receives contributions only from tree-level amplitudes and is 
suppressed by two factors: by the factor ~ (T w - W + /T lt z,Box) 2 from Eq. (|2.5|) , and by the 
additional factor ~ (m 2 . / m 2 Ni ) 2 comming from the heavy neutrino propagators. In Fig. 2, 
the choice = m^ 2 /3 was made since Majorana-type decays vanish if the masses of 
heavy neutrinos are equal. 

In the following, only the first group of decays is discussed. The results are given 
in Figs. 3-6. Figures 3 and 4 show the dependence of the branching ratios B(t~ — > 
e~ir + 7T~ /e~K + K~ /e~K°K°) on new weak interaction parameters of the model, and 
mjv^. Figures 5 and 6 illustrate the dependence of these branching ratios on model assump- 
tions for hadronic part of the amplitude and on some strong interaction parameters. 

The Figures 3(a) and 3(b) illustrate = = dependence of the branching 
ratios for (s v £) 2 = 0.01 and two different values of (s v £) 2 . The maximum values for branch- 
ing ratios are obtained for maximal m^, (s£ e ) 2 and {s v £) 2 values permitted by Eqs. (|3.1|) 
and 

B(t~ — > e~7r^7r~) ~ 0.74 • 10~ 6 (0.35 • 10~ 6 ) 
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B(T-^e~K + K-) ~ 0.42- 10~ 6 (0.20- 1(T 6 ) 
B(t~ -> e'K°K°) ~ 0.26- 1(T 6 (0.12- 10~ 6 ) (3.5) 

The expressions in the parentheses are obtained for the upper bound (s v £) 2 and (s£ T ) 2 
values refered in Eq. ( |1.4|) . The present experimental bound exists only for one of these 
decays 

B(t~ -> e-n + n-) < 4.4 • 10~ 6 , (3.6) 

because the main r~ — ► contribution mode to the r~ — > e~ K + K~ / e~ K° K° decays, 
r~ — > e~0, has not been experimentaly searched for yet. In Figures 3(a) and 3(b), the 
branching fractions B(t~ — > e~7r + 7r~ /e~ K + K~ /e~ K°K°) are shown. The behaviour of 
the branching ratio terms quadratic and quartic in s v £ expansion have similar behaviour 
as the corresponding terms in r — > e~M° decays F° r m N values below 200 GeV, 



quadratic (s^) 2 terms, that have ln(m 2 M / 'm^) large-m^v behaviour, prevail, while for larger 
rriN quartic terms having m 2 N large-mAr asymptotics dominate. As (s^ r ) 2 decreases, the 
branching fractions also decrease, but at the same time the pertubative unitarity upper 
bound onmjv increases, and, therefore, branching ratios increase in the larger interval. 
These two opposite effects lead to the small difference of the largest values for branching 
fractions in Eq. (|3.5| ). The nondecoupling behaviour of the branching ratios displayed 
in Figs. 3(a) and 3(b) is a consequence of the implicit assumption that the mixings s v £ 
may be kept constant in the whole m^-interval of interest. As mentioned in Introduction, 
s u £ oc m D /m M oc m D /m N%1 and, therefore, the constancy of s v £ implies that for large m N . 
values, the Dirac components, mo, are large also. Since the Dirac-mass values are bounded 
by the typical SM SU{2) x U(l) breaking scale, v ~ 250 GeV (more precisely, pertubative 
unitarity upper bound on the Dirac mass is < 1 TeV [[37]), this condition cannot 
be satisfied in the — ► oo limit, leading to vanishing effects of heavy neutrinos EC 



Nevertheless, for 0.1 TeV< < 10 TeV it can be fullfilled. Nondecoupling effects of the 



heavy neutrinos were first studied in Ref. []14| , and were also extensively studied in Refs. 

Figures 3(c) and 3(d) present the dependence of the branching ratios on {s u £) 2 and 
(s u £) 2 respectively, for = 4000 GeV. The branching ratios are almost quadratic func- 
tions of (sx, T ) 2 , and almost linear functions of (s£ e ) 2 . Such dependence is expected from the 
large-m^v behaviour of form factors [0] (see also Appendix B). 

Figure 4 illustrates Majorana-neutrino quantum effects. It displays the dependence of 
branching fractions on the ratio ttin 2 /m^ for fixed values = 1 TeV and w,n x — 0.5 TeV. 
The maximal B(t~ — > e~ -k+h^ / e~ K + K~ / e~ K° K°) values are obtained for m^/m^ ~ 3. 
These effects are also a consequence of large s u £ mixings (large Dirac components of the 
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neutrino mass matrix), since they enter through the loop functions depending on two heavy 
neutrino masses, which can be found only in quartic terms in the expansion. A similar 
behaviour has been found for t~ — > l'~M° and t~ — > I'^l^l^ jT2| decays. 

Figures 5(a)-5(c) show the influence of the main ingredients of the hadronic part of 
the amplitudes discussed in the comments of Section 2 on the branching ratios. Thick 
lines in Figs. 5(a)-5(c) correspond to the situation when one of the theoretical assumptions 
is changed. Thin lines serve as reference results and they coincide with the complete- 
calculation graphs shown in Fig. 3(a). 

Figure 5(a) shows the dependence of the branching ratios on vector- meson resonances. 
When the vector-meson propagators are replaced by their zero-momentum-transfer values, 
that is when the normalized vector-meson propagators (|2.27|) are replaced by 1, one obtains 
the chiral-limit values for the branching ratios plotted in Fig. 5(a), which are considerably 
smaller. The B{r~ — > e~7r + 7r~ /e - K~ K + ) branching ratios decrease by factors ~ 5 and ~ 
20 respectively. The decrease of r~ — > e~K~K + branching ratio is more prominent, because 
it receives main contribution from the narrower <p resonance, while to B(t~ — > e~7r + 7r~) 
only p-resonance contributes. The r _ — > e~K°K° branching ratio becomes almost equal to 
zero because its amplitude is proportional to the expression F^'Jf — which is almost 

equal to zero. 

In Figure 5(b), the U(3) x U(3)n/U(3)y breaking effects are emphasized by compar- 
ing the branching ratios obtained in the U(3) x U(3)n/U(3)y symmetric chiral model with 
reference results which include U(3) x U(3)n/U(3)y symmetry breakings. The symme- 
try breaking does not influence B(t~ — > e~7r~7r + ), but B(t~ — > e~ K~ K + /e~ K°K°) are 
enlarged by a factor ~ 1.5. 

The reference results include the #y- va hie derived from the Gell-Mann-Okubo 
quadratic mass formula, Qy = 39.1°. In Fig. 5(c), these results are compared with branch- 
ing ratios evaluated for 9y = 30°. As 9y is known to be close to the ideal mixing value 
arctan(l/\/2), the weak 6V-dependence displayed in Fig. 5(c) implies that 9y variation 
cannot influence the branching ratios strongly. 

The influence of the replacement ( |2.29|) induces so small changes of the branching 
ratios that they cannot be observed in a figure. For that reason these results have not been 
plotted. 

The Figure 6 gives the dependence of the partial decay rates on the momentum 
transfer variable, t = {p—p') 2 - The r~ — > e~7r~7r + decay rate receives the contribution from 
the broad p°-resonance only. The r~ — > e~ K~ K + je~ K°K° decays receive contributions 
from all three flavour-neutral resonances, but for the kinematical reasons only very narrow 
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^-resonance can be noticed in the spectra. 



4 Conclusions 

This paper completes the analysis of the experimentally investigated neutrinoless r- 
lepton decays within heavy-Majorana/Dirac-neutrino extensions of the SM, started in 
the previous publications [|l^,[n|. For the experimentally most promising decays, t~ — > 
ern + K~ j jjTTr" K + / fi + 7t~ K~ , the calculated branching ratios were found to be much 
smaller than the current experimental upper bounds. Nevertheless, the three of seven- 
teen explored decays, t~ — ► e~7r + 7r~ /e~K + K~ /e~K°K°, were found to have branching 
fractions of the order of 10~ 6 , and the first of them the branching fraction close to the 
current experimental sensitivity. The other two decays have not been measured yet, be- 
cause the reaction r~ — > e~0, giving the main contribution to these decays, has not been 
experimentally investigated yet. 

The main feature of the leptonic sector of the model used here is largeness of the 
heavy-light neutrino mixings ■ From it the dominance of the quartic terms and the 
m 2 N . behaviour of r~ — ► e~7r + 7r~ /e~K + K~ /e~K°K° in the large-mjy limit follows, giving 
rise to the enhancement of the branching ratios by the factor 40 relative to the results 
obtained by the analysis in which the respective terms are omitted. The behaviour 
and the m N , 2 /m Nl dependence of the branching ratios are also consequences of large s 1 ^ 
mixings. Particularly, the m^ 2 /m^ 1 dependence leads to the maxima of branching ratios 



for m^/mNj ~ 3, the same as in r — ► V l x l\ [12] and r — ► V M [ 13 1 decays. 



Several ingredients of the hadronic part of the r~ — > V~PxPi amplitudes, that influ- 
ence the magnitude of the corresponding branching ratios, were discussed. The most promi- 
nent contribution comes from the vector-meson resonances, giving rise to enhancemens of 
B{t~ — > e~7r + 7T~ /e~K + K~) by factors ~ 5 and ~ 20 and making B(t~ — > e~K°K°) 
different from its chiral limit value, zero, and approximately equal to branching values 
of the other two decays. The narrower resonances lead to larger enhancements. The 
U(3)ixU {3)r/U(3)v breaking of the chiral symmetry induce smaller changes of the branch- 
ing ratios, and they influence only the r~ — > e~K + K~ /e~K°K° branching fractions. All 
other modifications of or changes in the hadronic part of the r _ — > / /=F P 1 P 2 amplitudes 
discussed here have negligible influence on the branching ratios. 
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A Strong interaction Lagrangians 

The gauged chiral U{3)l x U(3)r/U(3)v Lagrangian extended by hidden U(3)i oca i 
symmetry and the mass term for pseudoscalar mesons reads 

C'A d£,y ~\- £ m ass ^kin 

= --jlTr{D^ L - D^ R ) 2 - jfiTriDptel + D^&f + + Am 

£2 , f2 



{^Tr{d^Jd^)) + U^rTr{m{U + 
+ {-e(a^ 2 )(p2 + ^ + ^)^ 

v yJnJ \ 2s w c w Hfl Vv/3 2s w c w ^2s w c w > 

rSy\- 2s 2 w C V 1 \ 
+ 7=- 



2swCw \/Q2swCw 

+V3s v uj^K + d a K- + K°d^K°) + VScv^K+d^K- + K°duK°) 

+K°*'»( - V2n + d^K- + n° du K° + V3c P K°duV + V^spK ^)^') 

+A°*' M (V2tt-^ K + - 7r°i K° - V3c P K°d^ V - V3s P K°d a //)} + ... (A. 1) 

where Chin is the kinetic Lagrangian of gauge fields, f n is the pseudoscalar decay constant, 
a is a free parameter equal 2 if vector-meson dominance is satisfied, g is the coupling of 
(hidden symmetry induced) vector mesons V, to the chiral fields £,l,r, cw — cos9\y, 

DM*) = (^-iV^(x))a(x)+i^(x)£ M (x), (L^R, C.^TZ,), (A.2) 
£ L>R ( X ) = e M-)//. e TM*)// % <j(x) = Q, (A.3) 
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a(x) = beeing special (unitary) gauge choice. £ M I 
gauge fields, 



x 



and TZ^x) 



are combinations of 



Ti^x) 



eQ[A^(x) - twZ^x)) + 



iQiA^x) - t w Z^(x] 



SwCw 



T z Z,(x) + -t=—W, 
V 2sw 



i-i- 



(A.4) 



where 





{ 2 ^ 


1 


f 1 \ 


Q = I 
V 3 


0-10 


' * ~ 2 


0-10 




Vo o -ij 




V -1 ) 



(A.5) 



are quark charge and isospin matrices, 

/ 



W^x) 



w, 



W-(x)c c 



,,, ,x)c c W+{x)s c \ 




V w- 



[X S r 











(A.6) 



c c and s c being cosine and sine of the Cabbibo angle, respectively. A^(x), Z^{x) and W^(x) 
are photon, Z-boson and W^-boson fields. The dots in Eq. ( |A.1| ) represent remaining terms 
in the gauged chiral U{3)l x U{3)r/U{3)v Lagrangian containing hidden U{3)i oca i symme- 
try, not interesting for the topics discussed in this paper. The first curly bracket contains 
minimal non-gauged chiral model Lagrangian. Using the Gell-Mann-Levy procedure [| 
the pseudoscalar-meson vector currents may be derived from that Lagrangian: 



V;(x) = -2Tr{T>(x),«Vr(x)]}, 



(A.7) 



with 7r(x) and T a defined below Eq. ( |2.16| ). For instance the vector current having quantum 
numbers of p meson reads 



-> 3 

V2 " 



1 



1 



1 



(A.8) 



V2 2V2 2V2 

Pseudoscalar mass terms may be found in the second curly bracket. The m is a mass matrix 
of u, d and s quarks, and r is defined in Eq. (|2 .17 ). Terms in the third curly bracket repre- 
sent photon-vector-boson and Z-boson-vector-boson interactions. These interactions define 
the corresponding gauge-boson-vector-meson coupling strengths (for instance photon-p- 
meson couping is equal to —eagf%). The fourth curly bracket comprises vector-meson-two- 
pseudoscalar-meson interactions and defines the corresponding couplings. 

The breaking of the U(3)l x U(3)r/U(3)v symmetry is introduced in the way of 
Bando, Kugo and Yamawaki [BO]. Besides the terms containing only the ^ or £r fields, 
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they added the additional mixing terms, combined with the matrix-valued parameters, 

( ° ^ 

e Ay = (A.9) 

V C A y J 

defining the magnitude of the symmetry breaking. These additional terms change the 
kinetic part of the pseudoscalar-field Lagrangian. To restore the original form of kinetic 
terms pseudoscalar-meson fields have to be renormalized: 

n(x) -> n r (x) = \ A /2 ir(x)\ A 2 (A.10) 

where \ A ,v = 1 + £a,v- Following the described procedure, one obtains the following 
expression 

pbr _ rbr , pbr , r> , r 

*~*A ' UlJ ~ , V ' *-mass ~r *-kin 

= '{- ±ffTr((D£ L & + D^ L s A e R ) - + D^ R e A £} L )) 2 

- a -flTr{{D,i L i{ + Dptev&t) + (D^r + D^ R e v ^ L )) 2 } + h.c. 
+C 

mass '-'kin 

2(1-7^) \ \ A , f , 2 ^ l-24, o 

2. ... l 



+ 1 

+ ! 



( Cv 1 U -2 _ 2s 2 1 + 7 M + St/ 1 (V - 2s 2 '- 



2 2(1-7^)- 



v^6 2swCw 



)k 



1^2^^ ZSw 3 ^ v^^cA 7 " 



— {' 



7T <9 M 7T~ 



7a 



3 ^ v/6 2svkCvk 



2(1-7^) 



)) 



V2V/3 



a + — -cv)(ir^x--x%x°) 

fo 1 + 2 7 2 7a 1 ) + + 2) - Cy{ lA l + 27V7A 1 )) 

-^((Tl 1 - lllA 1 ) + ^(lA 1 " 1) " Cvir/ - IVl'A 1 ))) 

x(K + duK'+K°daK°) 

+^[(^((7a 1 + 2 7 ^7l 1 ) + ^{l~A l + 2) - C v { lA l + 27V7A 1 )) 
+%^a - IvIa 1 ) + ^{Ia 1 " 1) " Cy{>y? - W))) 
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1 



1 

72' 



2v/3 



l + 2~/2 1 )K d f i(c P r ] + spr ) ') 



+ (--^)7 
a 



2 



(c P 77 + sp^d.K 



2 T 



v/6 



A ' spry + Cj V)0„Jr° ) + ^^^^^(cpTy + W) 



a i 

_3 _1 _3 1 

IvIa 7A 2 x r >oo / „ , „ a , /7f7a 2 7a' 



xK°d^c P ri' - s P rj) 



)K d II {c P 7] + S P 7]') + 



VE y/E 



+K 



0*,/* 



1/ 1 



7.7i(- 7r -^-^ r 



l + 2 7A 1 )ir%(cpr ? + spV) 



-1°^° - -1(1 - Ta 1 )^ (-^ + CpV; 



+(— - CV) 7a 5 - ^^ + 5 M 7T- + -K°d, 



27; 



v/2 



-sprj + cp^K ) - 
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-1 



{cpr] + s P r]')d^K c 



K°d^c P T] + sprf) 



7v7a 7a' 



a/3 2\/3 



\K°dJc P rj + Sprj') + 



'IVlA lA 



xK^J-sprj + Cprj'] 



+ 



(All) 



where ja,v — Ca,v + 1 and it, rj,. . . are renormalized pseudoscalar fields (superscript r is 
omitted). In the above expression, only the gauge-boson-vector-meson (first curly bracket) 
and vector-meson-two-pseudoscalar-meson (second curly bracket) interactions are kept. 



B Form factors and loop functions 



The composite form factors for r — > l'V° decays, a M o, b M o, and, c M o, appearing in 
the first three Eqs. ( j2.21| ), may be decomposed into the composite loop form factors F^ 1 ' , 



F T B l ' d x adb and F^ u in the following way 



n Tl ' 

a v o 



rriyo 



16M2 



w 7v° 



,Z rprl' 1 „Box,uu prV 



ayoFz +a 



I/O 



Box 



uu 1 „ Box,dd jTiTl'dd 



Box 



+a 



Box,ss 
V° 



rprl'ss , 
r Box a V° 



Box,ds rprl'ds , „ 



Box,sd jprl'sd 
r Box 
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h rl> 

0\/o 



4'o 



ia w m v o v 



16M^ 7y0 • 



(B.l) 



The factors «yo, /3yo and 7 y o, containing information on quark content of a vector meson 
V° (see Table II), and in part information on quark-7 and quark-Z° couplings, may be 
found in Table I. 

The loop form factors FT 1 ', G? 1 ' , Ff\ F Bo d x adb and F^ u , and Fj}' and G H ' contain 
the leptonic part of T 7 , Tz, Tbox and T# amplitudes, and may be further decomposed into 
elementary loop functions F 1 , G 1 , Fz, Gz, Hz, Fbox, Hbox, Fh, Gh, and Hh- The loop 
form factors FT 1 ' , GZj' , FT, 1 ', F Bo dadb and Fp^ u together with the elementary loop functions 



F 7 , G 7 , Fz, Gz, H z , F Box , H Box may be found in Refs. |12| , |T3|| . The composite loop form 
factor G H and the loop functions Fh and Gh were calculated for case of degenerate heavy 
neutrino masses in Ref. |14|]. Here the expressions for the composite loop formfactors Fj}' 



and G H are listed 



t?tV 



B* Ti B V j SijF H {Xi) + C* Gjj(Ai, Xj) + CijH H (Xi, Xj) 



I'Ni 



NiNj 



S NiNj (F H {X Ni ) ~ F H (0) + G H {\ Ni ,0) + G H (0, X Ni ] 



''rC'x N [GH{X.Ni, Aiv ) — Gni^Ni, 0) — Gjj(0, Ajv )) + C^N^HniXNi, X 



E F>* Ti Biij SijF H (Xi) + C*jGn(Xj, A») + CijH H (Xj, Aj) 



]T S^S^. $NiNj (F H (X Ni ) - F H (0) + G H (X Ni ,0) + G H (0, X Ni ] 
+Ch. n (Gif(Ajv,-, AjvJ — Gh(Xn p 0) — £?#((), AjvJ J + CN l N j H H {XM j , X^ i 



(B.2) 



together with the loop form factors G//, and contained in them 



G H {x,y) 



1 — x + xlnx/x^ xXh 



3 xlnx\x 

'2 + 1^)2 



3 xA# 

2 r~ 



(1 -x) 2 

1 — 4x + 3a; 2 — 2x 2 In x 

+ 2(1 -x) 3 

x(x — y)(l — x)(l — y) + x(l — y)(x + xy — 2y) lnx + y 2 (l — x) 2 lny 

-2(l-x) 2 (l-?/)(x-y) 2 

x In x — y In y — xy (In x — In y) 
x (x + y + xy) + — r y 

(1 -x)(l -y){x-y) 
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3 | (1 + x) lnx — (1 + y) lny ^ 1 / In a; \ny \\ 



4 2(x-y) 2(x-y)^-l+x -1+y 

TT , , , f x(x — y)(l — x)(l — y) + x(l — y)(x + xy — 2y) lnx + y 2 (l — x) 2 lny 

////(•'••//! = v 7 ^ " 



V -2(l-x) 2 (l-y)(x-?/) 

,„ 1. xlnx — ylny — xw(lnx — lnw) 

x(2 + -(x + y)) + 



2 



2 V 

3 (1 + x) lnx — (1 + y) lay 1 / In a; \ny 

4 + 2(x - y) + 2(x -y)\-l+x -1 + y 



(B-3) 



For reader's convenience, Fh, Gh, and if// are evaluated for some special values of argu- 
ments, 

Fh(0) = ~ F H {1) = ^, 

— 5x + Ax 2 + x 3 — (10x 2 — 6x 3 + 2x 4 ) lnx 
G ^ (X ' X) = 4(1^ ' 

-3 + 17s - 13a; 2 - x 3 + (14a; 2 - 2a; 3 ) In a; 



G H {x, 1) 
G H (l,x) 



4(1 -a;) 3 

1 - 7x + 8x 2 - 5x 3 + 3x 4 - (6a; 2 - 2x 3 + 2a; 4 ) In x 



4(1 - x\ 



3 



„ , — x + x 2 — xlnx _ , — 3x + 2xlnx 

G H {x,0) = 2(1 -x) 2 ' G ^(°' a; ) = 4 . 

= G H (0,0) = 0, G h (0,1) = ~ G? ff (l,0) = i 

rr . . -5x + Ax 2 + x 3 - (10x 2 - 6x 3 + 2x 4 ) lnx 

H H (x,x) = 



4(1 -x) 3 
xi(7 - 8x + x 2 + (3 + 4x - x 2 )lnx) 
2(1 -x) 3 ' 
£2 (— 5 + 7x — llx 2 + 9x 3 — (8x — 2x 2 + 6x 3 ) lnx) 



H H (x,l) 

Hh{1 > X) 8(1 -x) 3 

H H (0,x) = H H (x,0)=H H (l,l) = 0. (B.4) 

If s'i are kept constant, all composite loop form factors are increasing functions of the 
heavy neutrino masses. The asymptotic behaviour of the form factors F^ 1 , G'L 1 and F^ 1 ' , 
in the limit Ai ^> 1 and p — A2/A1 > 1, are listed in Ref. fll2| . Here we list the form factors 
F]j' and in the same limit, 



+P 1 ' 2 

tig 



+ s^' Vfs^ 2 3pA 1 (4 + 4pV 2 + (l-p 1 / 2 )lnp) 

+SiSi y Si) 4(1 + P V2)3 • ^ 
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C Phase space functions 



The momentum dependent part of the absolute squares of the t~ — ► Z /=F PiP 2 am- 
plitudes may be expressed in terms of the Mandelstam variables t = (p — p') 2 and 
s i = ip' + Pi) 2 = (p ~ P2) 2 ■ The r~ — > l' T PiP 2 decay rates contain the integrals of 
the corresponding absolute squares of the amplitudes over si and t variables: 



(m—m 1 ) 2 



dt f 1 ds 1 {\T{ T -^l'*P 1 P 2 )\ 2 ) (CI) 



2567r 3 m 3 J( mi +m2) 2 

where (|T(t _ — > l' T PiP 2 )| 2 ) is the square of the amplitude averaged over initial and summed 
over final lepton spins. The boundary si-values, sf(t), read 



± M 2 2 B{t) ^^B(tf-AA(t)C(t) 
s f(t) = m 2 + m 2 + -r4-f ± v 



A{t) 



A{t) 



(C.2) 



where 



A(t) 
Cit) 



it, B(t) = -2(m 2 -m' 2 + t)(t + m 2 2 -m\) 
m 2 (t + m\ — m\) 2 + m\\{m 2 , m' 2 , t), 



(C.3) 



and X(x, y, z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz. Since the momentum dependent parts of 
the squared amplitude in Eq. ( |(J.1| ) contain only powers of the Si variable, Si integration 
is easily performed resulting with expressions which are denoted by a, (3, 7, 5, e, (, 77, 
1, K and ui: 



a 

7 
5 



2S\ + S\ 



2t - 2(m 2 + m' 2 + m{ + m 2 2 ) 
1 



b 1 2 1 /2 \ 1 / 2 1 /2\2 1 o 2 2 

-(m +m ) + -{m + m ) + 2m l m 2 



Si 

s? 
si 



2 12 

m — m 



t I 2 i2\ 1 / 4 /4\ / 2 2 /2 2\ 

-(m — m J (m — m J — (m m 1 — m m 2 ) 



t 1 2 1 /2 \ 1 / 2 /2\2 

-(m +m J + — (m — m ) 



/ 2 I2\i 2 2\ 

(m — m )(m 1 — m 2 J 



+ 5? 



^ / 2 1 /2\ 1 l/„2 /2\ / 2 2\ 

— -(m +m J + — (m — m )im 1 — m 2 ) 



1 1 / 2 /2\2 1 / 2 1 2\/ 2 1 /2% 1 1 / 1 / 4 /4\/ 2 2\ 

+-(m — m ) + {m 1 + m 2 ){m +m J + ( — ~ m )\ m i~ m 2) 

I 2 2 12 2\i 2 2\ ( 2 1 2\/ 2 /2\2\ 

— imm x — m m 2 )\m x — m 2 ) — [m 1 + m 2j )(m — m J I 



5? 



-(m + m ) 



/ 4 /4 \ / z z\ 

- — (m — m — m 2 j 



2(m z + m") - 
f 



'2) _ 2 A 2 + m /2 )2 + (m 2 + m '2 )(m 2 + m 2^ 



2 1 /2 \ / 2 12 \ 2 

-(m +m ) — — [m — m ) 



27 



n 



i - 

K = 
UJ = 

where 



— (m 2 + m' 2 ){m\ + mj) + -(^2m 2 m' 2 (m 2 + m' 2 ) — Am 2 m' 2 {m\ + m 2 ) 

I 4 IA\/ 2 2\ , ^ / 2 , /2\/ 2 , 2\2\ 

— (m — m J(m 1 — m 2 ) + -(m +m ){m 1 +m 2 ) j 

, 1 /o/ 4 /4\/ 2 2 /2 2\ , / 4 /4\/ 4 4\ 
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Si 



2t 



/„ 2 /2\/ 2 , /2, 2, 2\ 

(m — m )(m +m + m 1 + m 2 j 



i / 2 /2 i 2 2 \ 

+-(m — m +m 2 — m 1 J 



2 

mm bi 



m 2 + m' 2 — t) 



s? 



f i 2 2\2i 2 , /2 j.\ 

-{t — m x — m 2 ) [m + m —t) 



(C.4) 



(C.5) 



The definitions of other quantities in Eqs. (|C.1| - |C.3|) may be found in the previous text. 
The t-integration of expressions (|C.1| ) has been performed numerically. 
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Figure Captions 



Fig. 1: Feynman graphs pertinent to the semileptonic lepton-number- violating de- 
cays t~ — ► l' + P{P<2 [Fig.l(a)] and to the semileptonic lepton-flavour vio- 
lating decays r~ — > l'~PiP 2 [Fig. 1(b)]. The hatched blobs represent sets 
of lowest order diagrams contributing to three-point and four-point func- 
tions violating lepton flavour. These sets of diagrams may be found in Refs. 
I^ , |l3||L4] , pl^|l6fl . The double hatched blobs represent interactions through 



which the final state pseudoscalar mesons are formed. 

Fig. 2: Branching ratios (BR-s) versus heavy-neutrino mass = = \mN 2 for 
the decays r~ — > e~7T~7T + (thick solid line), r~ — > e~K~K + (thick dashed 
line), t~ -> e-K°K° (thick dot-dashed line), t~ -> e~ "7r~ K + / 'e~ "7r + K~ (1), 
T~ -> e-n°K°/e-iT K (2), r~ -> e~7]K° /e~r]K° (3), r" -> e~ rf K° / e~ rf K° 
(4), r~ — > e _ 7r°7r° (5), r~ — > e - ?]^ (6), r~ — > e~T]rf (7), r~ — > e + 7r~7r~ 
(8), r~ -> e+vr-fsT- (9), r" -> e + K~K~ (10), assuming (s£) 2 = 0.01 and 
( s ^)2 = 0.05. 

Fig. 3: Branching ratios versus new electroweak parameters of the model. Fig. 3(a): 
BR-s versus = = mjv 2 , assuming (s^ ) 2 = 0.01 and (s£ T ) 2 = 0.05. 
Fig. 3(b): BR-s versus = = rriN 2 , assuming (s v £) 2 = 0.01 and 

(s£-) 2 = 0.02. Fig. 3(c): BR-s versus (s^) 2 , assuming m N = 4000 Gel/ and 
(s u L e ) 2 = 0.01. Fig. 3(d): BR-s versus (s u L e ) 2 , assuming m N = 4000 GeV and 
(#) 2 = 0.05. 

Fig. 4: Branching ratios versus ratio mN 2 / m N 1 for the decays of Fig. 3, assuming 
m7Vi = mjV2 = 4 TeV, (s£ e ) 2 = 0.01 and (s u £) 2 = 0.05. 

Fig. 5: Branching ratios versus — = for the decays of Fig. 3, assuming 
( s l ) 2 = 0-01 an d ( s zT) 2 = 0.05. The figure illustrates the dependence of 
BR-s on few ingredients of hadronic part of the amplitudes. Fig. 5(a): The 
influence of the vector meson propagators on BR-s. Fig. 5(b): The influence 
of the £/(3)l x U(3)n/U(3)v breaking on BR-s. Fig. 5(c): BR-s for 9 V = 30°. 
Thin lines represent the reference graphs and coincide with thick lines in the 
Fig. 3(a). Thick lines show BR-s in a situation when one of the ingredients 
of the hadronic part of the amplitudes is changed. 

Fig. 6: Partial decay rates devided by the r decay width as functions of t — {p — p') 2 
assuming m Nl = m N2 = 3700 GeV, {s v £) 2 = 0.01 and {s u L T ) 2 = 0.05. 
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Table I: Coefficients defining composite form factors for r — > l'V° decays: 

Besides the constants listed in the Table I, there are two more constants different from 
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Table II: Quark content of the pseudoscalar meson states and fields: 

The meson states listed in the Table II correspond to the tensor description of meson states, 

what is more appropriate for chiral model calculations. The states \tt + ) and K°) have op- 
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